External magnetic field effects on a distorted kagome antiferromagnet 
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We report bulk magnetization, and elastic and inelastic neutron scattering measurements un- 
der an external magnetic field, H, on the weakly coupled distorted kagome system, Cu2(OD)3Cl. 
Our results show that the ordered state below 6.7 K is a canted antiferromagnet and consists of 
large antiferromagnetic ac-components and smaller ferromagnetic 6-components. By first-principle 
calculations and linear spin wave analysis, we present a simple spin hamiltonian with non-uniform 
nearest neighbor exchange interactions resulting in a system of coupled spin trimers with a single-ion 
anisotropy that can qualitatively reproduce the spin dynamics of Cu2(OD)3Cl. 



PACS numbers: 75.10.Jm, 75.25. +z, 75.50.Ec 

The quest for novel quantum spin states in a space 
higher than one-dimension has been active [l[ since the 
resonating valence bond state was proposed as the ground 
state of a two-dimensional quantum spin triangular sys- 
tem and was also proposed to be relevant to high T c 
superconductivity in cuprates Recently, a few com- 
pounds have attracted lots of attention as good model 
systems [jj . One of them is Cu4- x Zn x (OD)6Cl2 which 
was initially believed to realize a two-dimensional kagome 
lattice with quantum spins when x = 1 and a three- 
dimensional pyrochlore when x = 0.[5] The prospect 
looked promising because Cu4_ x Zn x (OD) 6 Cl2 does not 
order down to 50 mK for most values of x. Q 

Recent neutron diffraction studies, however, have 
found that Cu3Zn(OD) 6 Cl2 with rhombohedral i?3m 
symmetry realizes a perfect kagome lattice but with 10% 
nonmagnetic site disorder, [jj [8| Furthermore, when Zn 
ions are replaced by Cu ions, Cu4_ x Zn x (OD)gCl2 with 
x < 0.3, becomes monoclinic (P2i/n) due to the Jahn- 
Teller effects of Cu 2 +(3d 9 ) ion. This leads to different 
bond angles for different Cu-O-Cu superexchange paths 
in such a way that the kagome lattice is distorted to 
have nonuniform nearest neighbor interactions and the 
interplane coupling between the kagome and triangular 
planes is expected to be weak. Thus Cu4_ x Zn x (OD) 6 Cl 2 
realizes a good two dimensional quantum spin system for 
all x but it does have defects with regards to an ideal 
kagome quantum antiferromagnet: for x — 1 the non- 
magnetic site disorder while for x — the nonuniform 
coupling constants. 0] 

The distorted kagome quantum spin system 
Cu 2 (OD) 3 Cl is interesting in its own right. Spe- 
cific heat measurements have shown that upon cooling 
it undergoes two transitions, one at T C 2 = 18 K and 
another at T& — 6.7 K. [9( Muon spin relaxation (/iSR) 
measurements yielded evidence for long range order at 
T C 2 and for enhancement of strong spin fluctuations 



at T c i. [ij [l(| However, in the specific heat data the 
amount of entropy given off at the T C 2 transition is tiny 
in comparison to the change in entropy at T c \, which 
appears inconsistent with the /xSR results. Neutron 
diffraction studies revealed development of observable 
magnetic Bragg peaks only below T cl Q , which indicates 
that this transition is associated with Neel order, a 
result that is consistent with the large amount of entropy 
released at T c \ but not at T C 2- The spin structure was 
found to be collinear rather than the 120° configuration 
that is expected for a uniform kagome antiferromag- 
net. A recent theoretical study showed that collinear 
spin structures can indeed be favored by a distorted 
kagome. fi"lj| For T c \ < T < T c2 , on the other hand, there 
were no observable magnetic peaks but instead a weak 
magnetic excitation was observed at an energy transfer 
Klu = 7 meV.0] This intermediate phase was proposed 
to be a quantum spin singlet state, in contradiction with 
the /xSR results. 

Here we present our bulk magnetization and neutron 
scattering measurements on Cu2(OD)3Cl under an exter- 
nal magnetic field. Our principal results are the follow- 
ing. (1) The long-range magnetic order below 6.7K con- 
sists of strong antiferromagnetic collinear ac-components 
(| (M) ac |= 0.56(10) ^ B /Cu 2+ ) and a weak ferromag- 
netic 6-component (| (M) b |= 0.151 /i B /Cu 2 +). (2) The 
7 meV peak in the magnetic excitation spectrum, that 
was previously ascribed to a quantum singlet state, is due 
to a van Hove singularity of the spin wave excitations at 
the zone boundary of the excitation spectrum. (3) The 
spin dynamics can be qualitatively explained by a simple 
spin hamiltonian with non-uniform nearest neighbor ex- 
change interactions resulting in three different coupling 
constants with values of 8 meV, 3.7 meV and 2.8 meV in 
the kagome plane and a single ion anisotropy of 0.4 meV. 

A 20 g powder sample of Cu2(OD) 3 Cl was made at 
the University of Fukui using the standard method re- 
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FIG. 1: (a) Bulk magnetization, (M)bulh, and (b) suscepti- 
bility, d{M)buik/dH, measured as a function of an external 
magnetic field at several different temperatures. The inset 
shows the low field (-2 T < H < 2 T) data of (M) bu i k at five 
selected temperatures. 



ported elsewhere. [B[ Bulk magnetization measurements 
up to H = 12 T were performed at the Hahn-Meitner 
Institut (HMI), Berlin, while high field measurements up 
to H = 35 T were done at the IMR, Tohoku University. 
Neutron scattering measurements were performed at the 
NIST Center for Neutron Research (NCNR) and at the 
Berlin Neutron Scattering Center (BENSC) at HMI. De- 
tails of the neutron scattering experimental setups are 
described in the captions of the figures that present the 
data. 

Fig. 1 (a) shows bulk magnetization, (M)buik, mea- 
sured as a function of an external magnetic field, H , at 
various temperatures, spanning the two phase transitions 
at T cl = 6.7 K and T c2 ~ 20 K. At T = 1.5 K << T cl , 
(M)buik has a nonzero value of 0.15 /j,b/Cu 2+ at H — 
Tesla and exhibits hysteresis with H (see inset). This 
tells us that the Neel state below T cl has ferromagnetic 
components. As shown in Fig. 2 (a), below T c i, antifer- 
romagnetic Bragg reflections with the characteristic wave 
vector of Q m = (0,0,0) appear. The ferromagnetic com- 
ponent clearly observed in (M)b u lk would yield magnetic 
scattering on the top of the nuclear Bragg reflections 
only. We have applied the group theoretical approach 
in solving the magnetic structure. There are four irre- 
ducible representations of group P2\/n for Q m . Among 
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FIG. 2: (a) Elastic neutron scattering as a function of the 
scattering angle, 26. The spin structure projected (b) on the 
ac-plane and (c) on the afe-plane. Blue and green spheres 
in (b) and (c) represent kagome and triangular Cu 2+ spins, 
respectively. 



them, the representation t\ has antiferromagnetic basis 
functions for the crystallographic a- and c-components, 
ij)* = (1,0,0) and ip^ = (0,0,1), respectively, and fer- 
romagnetic fe-component, = (0, 1,0). We have found 
that the linear combination of the basis functions of the 
n representation, 0.55(5)^ + 0.151^ + 0.23(25)^, 
reproduces our elastic neutron scattering pattern (Fig. 2 
(a)). As shown in Fig. 2 (b), the ac-components form 
the collinear, antiferromagnetic spin arrangement that 
was proposed previously. [7J 

When an external magnetic field is applied up to about 
10 Tesla, {M)b u ik increases rapidly to the average value 
of 0.7 ^b/Cu 2 +. For H > 10 Tesla, {M) bulk keeps in- 
creasing at a slower rate to 0.9 /is/Cu 2+ at H = 35 
Tesla which is close to the value of 1 /is/Cu 2+ , ex- 
pected when all spins are fully polarized ferromagnet- 
ically. The change in the increase rate can also be 
seen in d(M)f Ju ik/dH vs H (Fig. 1 (b)) which exhibits 
a broad peak at around H — 8 Tesla and decreases 
rapidly at higher H, indicating a field-induced mag- 
netic phase transition. Upon warming, the antiferromag- 
netic ac-components weakenQ, and so does the ferro- 
magnetic 6-component (see the 6 K data in the inset), 
and both components disappear at T c \ = 6.7 k[|| 0]. 
At T = 6 K ~ T c i, the broad peak which was found in 
d{M)buik/dH at H c i ~ 8 Tesla at T « T c \ becomes a 
broader bump at a lower H c i ~ 5 Tesla. While in the in- 
termediate phase, T c i < T < T C 2, the broad peak found 
in d(M)b u ik/dH for T < T c i is now replaced by a sharp 
peak at low H; H c2 = 1.5, 2, 2.2 Tesla for T = 8, 10, 
12 K, respectively, and at 16 K a broad peak at around 
H c2 = 2 Tesla. At 30 K and 40 K > T c2 , d(M) bulk /dH is 
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FIG. 3: (a) Elastic neutron scattering pattern as a function 
of Q up to 1.3 A" 1 , obtained at T = 0.1 K under an external 
magnetic field of T (blue) and 3.5 T (red), (b) The differ- 
ence between the two diffraction pattern, I (H = 11 T) — I (H 
= T). The data were obtained with an incident neutron en- 
ergy of Ei — 3.7 meV and collimation of guide-80-80-open at 
the cold-neutron triple-axis spectrometer, SPINS, located at 
the NCNR. Contamination from higher-order neutrons were 
eliminated by placing a BeO filter after the sample. Inset of 
(a): H-dependence of the (010) Bragg peak intensity at four 
different temperatures, T = 0.12 K, 2.2 K, 4.3 K and 4.9 K. 
The measurements were done at the cold neutron triple-axis 
spectrometer V2 at HMI with Ei = 5 meV. 



nearly flat, i.e., (M)^^ increases linearly with H, up to 
H = 14 Tesla, suggesting paramagnetic behaviour. 

To study the field effects on the Neel state below T c \, 
we have performed elastic neutron scattering measure- 
ments with an external magnetic field. Fig. 3 (a) shows 
the elastic neutron scattering pattern as a function of mo- 
mentum transfer, Q, measured at T = 0.12 K with two 
different fields, H = T (blue symbols) and 11 T (red 
symbols). The two antiferromagnetic Bragg peaks at the 
(001) and (010) reflections, that are prominent in zero 
field, become very weak at H = 11 T, indicating a field in- 
duced phase transition occurred below H = 11 T. The in- 
set of Fig. 3 (a) shows that at T = 0.12 K the peak starts 
weakening at around 7 T which is close to the H c \ = 8 
T observed in d{M)buik/dH and disappears at around 
10 T. When temperature increases, the critical field de- 
creases, for instance H c \ = 6.5 T at T — 4.9 K, which 
is consistent with our bulk magnetization data. Where 
did the antiferromagnetic scattering go? The intensity 
difference between H = and 11 T, plotted in Fig. 3 
(b), shows that the peak at Q = 1.16A -1 which corre- 
sponds to the nuclear (1,0,-1) & (0,1,1) reflections has 
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FIG. 4: Inelastic neutron scattering intensity obtained at Q — 
0.6 A -1 under external magnetic fields of, H = T (blue), 5 
T (green), 9 T (pink) and 13.5 T (red), at temperatures of 
T = 1.5 K (a) and 10 K (b). The data were obtained at the 
cold-neutron triple-axis spectrometer V2. The data shown in 
the inset of (a) were obtained at Q = 1.2 A -1 , at the SPINS. 
At both spectrometers, the energy of the scattered neutrons 
was fixed to be Ef = 3.7 meV and a horizontally focusing 
analyzer was utilized to increase the sensitivity. 



increased upon application of the field. Altogether, the 
neutron and susceptibility data suggest a picture where 
below 8T the field increases the size of the spin moment 
without changing its direction significantly. In contrast, 
above 8T where the ferromagnetic component increases 
at the expense of the antiferromagnetic component and 
the spins are able to increase their alignment with the 
field perhaps by overcoming an anisotropy. 

Let us now turn to field effects on spin fluctuations. 
Fig. 4 (a) shows the H-dependence of the magnetic ex- 
citation spectrum measured at T = 1.5 K « T c \. At 
H = Tesla, there are two modes centered at TkJ = 1.32 
meV and 7 meV, as previously reported 0|. When H is 
applied, the sharp peak at 1.3 meV splits into two broad 
peaks and the splitting between the two peaks increases 
with H. At H = 13.5 T, the two peaks are centered at 
Tiui = 1.25 meV and 2 meV, respectively. The 7 meV peak 
also splits into two peaks (see inset), which unambigu- 
ously tells us that it is not due to the singlet-to-triplet 
excitations of a quantum dimer as previously thought, 
but is due to spin waves. At 10 K > T c \, and zero field, 
the Tiuj = 1.3 meV mode is replaced by a low energy con- 
tinuum (see Fig. 4 (b)), commonly found in a cooperative 
paramagnet. When the magnetic field is applied, distinct 
excitation modes appear out of the continuum. For H > 
9 T, there exist a broad feature that can be fitted by two 
gaussians. At H = 13.5 T, the two peaks are centered 
at 1.25 meV and 2 meV. This behavior is quite different 
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FIG. 5: (a) Dispersion of spin wave excitations along (h,0,0), 
(l,k,0), and (1,1,1) directions, (b) Powder- averaged neutron 
scattering intensity as a function of Tiuj. Black and blue lines 
are the results for H = OT and 11.5 T, respectively. The inset 
of (b) shows the kagome lattice with the ac-components of 
the moments. Red, blue, green lines represent different bond 
strengths. 



from the field-induced single-ionic Zeeman-likc behavior 
observed in Cu 4 _ x Zn x (OD) 6 Cl 2 (x=0.66 and 1)@,I3. In- 
stead, it is very similar to the behavior observed at 1.5 K 
(see Fig. 4 (a)) in an external field. This indicates that 
the paramagnetic continuum at low energies above T c \ is 
due to spin fluctuations anticipating the canted antifer- 
romagnetic order below T c \. This also suggests that the 
7 meV mode that survives above T cl might be related to 
the spin wave fluctuations of the Neel state. 

It is almost impossible to extract the exact spin hamil- 
tonian, TL 1 from powder-averaged inelastic neutron scat- 
tering data, because the crystallographic directional in- 
formation is lost. This is especially the case when 7i 
is complex as in Cu2(OD)3Cl. Therefore we have per- 
formed the local-density approximation (LDA)+[7 (the 
d — d Coulomb interaction parameter) calculations for a 
wide range of U to obtain the ratios of different exchange 
constants. Our results tell us that there are three differ- 
ent coupling constants in the kagome plane as shown in 
the inset of Fig. 5b with the ratios of J k / J k ~ 0.46 and 



J k I J k — 0.35. [12| The inter-plane couplings turned out 



to be much weaker, less than 0.1J fc . We then performed 
linear spin wave analysis with those ratios fixed and with 
two fitting parameters, J k , that couples the spins into 
trimers, and a single ion anisotropy along the spin di- 
rection, D, for the collinear AFM structure composed of 
the ac-components of the Cu 2+ ions, neglecting the FM 
^-components. 

We focused on the following experimental observations: 
upon warming the 1.3 meV mode shifts to lower energies 
as the frozen moment weakens while the 7 meV mode 
does not shift in energy. This suggests that the 1.3 
meV excitation is a zero energy mode shifted to finite 
energy by anisotropy, while the 7 meV mode is due to 
the van Hove singularity of the excitations at the zone 
boundary of the excitation spectrum. As shown in Fig. 
5 (b), a sharp and strong peak at 1.3 meV and a weak 



and broad peak at 7 meV can be generated in zero field 
with J k = 8 meV and D = —0.4 meV. Fig. 5 (a) show 
that the 1.3 meV peak is due to two nearly-flat doubly 
generate modes, while the 7 meV peak arises from the 
enhanced density of states at the zone-boundary energy 
of another doubly degenerate mode. The fourth doubly 
degenerate mode generates two broad peaks around 4 
meV, which have not been observed experimentally. The 
discrepancy may tell us that the exchange interactions 
are not isotropic or the asymmetric exchange interactions 
such as the Dzyaloshinskii-Moriya interactions need to be 
included in our analysis. 

When a Zeeman term, g^sH ■ ^2^?, is included in 
the spin hamiltonian, the 1.3 meV peak splits into two 
peaks at a lower and a higher energy, while the excita- 
tions around 7 meV becomes broader and exhibit two 
peaks at around 6 meV and 7 meV, which are qualita- 
tively consistent with the data shown in Fig. 4 (a). Thus, 
the previous interpretation of the 7 meV excitation as 
a triplet associated with spin dimcrs in the intermedi- 
ate temperature (6.7 K < T < 18 K) regime Q should 
be modified to reflect the large intratrimer interactions. 
So far, the observed characteristics of this phase are the 
weak anomaly in the specific heat, /iSR measurements, 
and the small kinks in our d(M)b u i k /dH vs H data. Fur- 
ther experimental and theoretical investigations are nec- 
essary to understand the nature of the intermediate T 
phase. 

In summary, using bulk susceptibility, and elastic and 
inelastic neutron scattering measurements with and with- 
out an external magnetic field, we have characterized 
the static and dynamic spin correlations in Cu2(OD)3Cl. 
Our analysis shows that Cu2(OD)3Cl should be consid- 
ered as a system of coupled spin trimers on the kagome 
lattice. 

We thank A. B. Harris, S. Nagler, T. J. Sato, Y. B. 
Kim, M. Kofu for helpful discussion. 
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